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We study pinball billiard dynamics in an equilateral triangular table. In such dynamics, collisions with the walls are non- 
elastic: the outgoing angle with the normal vector to the boundary is a uniform factor A < 1 smaller than the incoming 
angle. This leads to contraction in phase space for the discrete-time dynamics between consecutive collisions, and 
hence to attractors of zero Lebesgue measure, which are almost always fractal strange attractors with chaotic dynamics, 
due to the presence of an expansion mechanism. We study the structure of these strange attractors and their evolution 
as the contraction parameter A is varied. For A G (0, ^), we prove rigorously that the attractor has the structure of a 
Cantor set times an interval, whereas for larger values of A the billiard dynamics gives rise to nonaccessible regions 
in phase space. For A close to 1, the attractor splits into three transitive components, the basins of attraction of which 
have fractal basin boundaries. 



Billiard models form an important class of dynamical sys- 
tems, in which a particle collides with fixed boundary el- 
ements. The collision rule plays a crucial role in the dy- 
namics: for classical elastic collisions, phase-space volume 
is preserved, and the dynamics is typical of Hamiltonian 
systems^. The recently-introduced pinball billiards^, or 
pinbilliards^ instead have a non-elastic collision rule, remi- 
niscent of the effect of the extra impulse in the normal di- 
rection received by a ball in a pinball machine when it hits 
an active bumper. The present authors previously studied^ 
a rule in which the outgoing angle at a collision, with re- 
spect to the normal vector to the boundary, is shrunk by 
a factor A < 1 compared to the incoming angle. In con- 
vex billiard tables with focusing and fiat boundary com- 
ponents, attracting periodic orbits, strange attractors, and 
bifurcation phenomena are then observed, depending on 
the boundary geometry. When the curvature of a bound- 
ary element is non-zero, some hyperbolic (chaotic) prop- 
erties are obtained, both in the elastic^ and non-elastic^ ^ 
cases. 

In order to make further progress in understanding the 
origin of the complicated behaviour induced by the non- 
elastic collision rule, in this paper we restrict attention to 
pinball billiards in polygonal tables, specifically the equi- 
lateral triangle, which turns out to capture many of the 
main features of such systems. Here we again find hyper- 
bolic (chaotic) strange attractors, whose structure evolves 
in complicated ways as A varies between and 1. Due 
to the simplicity of the geometry, however, many of these 
structural changes may now fact be understood analyti- 
cally. 
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Polygonal billiard tables, in the classical case of elastic colli- 
sions, exhibit dynamical properties which depend strongly on 
the angles of the table: if the angles are rational multiples of 
TT, then the dynamics are completely regular, and only a few 
angles can be realised during a given trajectory. On the other 
hand, arbitrarily close to a given polygon is another (with ver- 
tices as close as desired), whose billiard map is ergodicP. Fur- 
thermore, triangular billiards whose angles are all irrational 
with TT are conjectured to be mixing^, although there is cur- 
rently no proof of this property^. 

The specific non-elastic rule that we stud> ^ shrinks the outgo- 
ing angle at a given collision by a uniform factor < A < 1 
with respect to the incoming angle at that collision, mea- 
sured from the normal vector to the boundary. This gives a 
one-parameter family of maps interpolating between a one- 
dimensional map for A = and elastic billiard dynamics for 
A = l. 

We might expect that the resulting phase-space contraction 
would lead only to asymptotic periodic orbits. In fact, how- 
ever, there is also a strong expansion mechanism present, and 
the interplay between these gives rise to hyperbolic (chaotic) 
strange attractors in phase space. 

This possibility was foreseen in a result of one of the present 
authors, together with Pujals and Sambarino^, that all pinbil- 
liards on polygonal tables have dominated splitting, which is 
an extension of hyperbolicity. Pujals and Sambarino^ gave 
a trichotomy for invariant sets with dominated splitting on a 
compact manifold: they must be either finitely many intervals 
of periodic orbits with bounded period; finitely many simple 
closed curve where the dynamics is conjugate to an irrational 
rotation on a circle; or a finite union of hyperbolic sets. 
In polygonal tables, the pinball billiard map is only piecewise 
continuous, and hence the phase space is non-compact. For 
the equilateral triangle, we find a non-compact attracting set 
for any A. Our objective is to study the structure and dynamics 
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FIG. 1. (a) Non-elastic reflection rule; (b) several bounces in an 
equilateral triangle for A =0.5, with collisions numbered in order 
starting from the initial position 0. 



of these attractors. The results of the present paper suggest 
that for any < A < 1, the attractors exhibit a rich fractal 
geometry. In fact, for A < ^ we rigorously prove the existence 
of a transitive attractor whose structure is roughly a Cantor set 
of lines, and give a symbolic model that completely describes 
the dynamics. Moreover, the attractors for all values of A in 
this regime are topologically conjugate to one another. 
For larger values of A, the attractors undergo a sequence of 
structural changes, giving rise to striking images, shown in 
fig. [2j Some regions of phase space become inaccessible due 
to the nature of the billiard dynamics, which is reflected in 
the geometry of the attractor. A hyperbolic periodic orbit then 
becomes isolated from the transitive strange attractor, and fi- 
nally, for A close to 1, the strange attractor splits into three 
distinct transitive components. 



I. DEFINITIONS AND MAIN RESULTS 

Let 7 be the boundary of a billiard table, taken to be a closed, 
connected, convex curve in two dimensions, with length |7|, 
and certain smoothness conditioners xhe boundary of the 
table is parametrised by arc length r. 

We consider non-elastic billiards in which the angle of reflec- 
tion is modified from an elastic one. General modifications 
were previously studied^; in this paper we restrict to the fol- 
lowing one. The pinball billiard (pinbilliard) map with 
parameter A is given by 7\(^07 ^o) = (n, ^i), where n is ob- 
tained as in the elastic billiard, by moving along the direction 
determined by the outgoing angle Oq, starting at the boundary 
point ro, until the next boundary collision. The new outgo- 
ing angle 6i (at the next collision) is then given, following a 
standard sign convention, by 0i := — At]!, where rfi is the an- 
gle between the incident velocity vector at qi and the outward 
normal —n{qi)', it satisfies — ^<Q\<^. This is depicted 
in figure [T] together with an example trajectory in an equilat- 
eral triangle. 

The phase space of the pinbilliard map is thus contained in 
M := [0, |7|) X (-|, |). In fact, there is a finite number of 
curves where the map is not defined, consisting of points 
{s^B) located at a vertex of the polygon, or whose trajecto- 
ries hit a vertex in one iteration, so that an outgoing angle 
cannot be uniquely defined. The singular locus Sing(A) C M 
of Tx consists of all preimages of these curves, for which 



the dynamics can no longer be continued after some iterate 
at which the dynamics hits a vertex. We define the phase 
space M := M\ Sing(A), the set of points which can be it- 
erated arbitrarily many times. Note that Sing(A) has zero 
Lebesgue measure, since it is the union of countably many 
regular curves. 

The main object of interest is the structure of the maximal for- 
ward invariant set in M, which we refer to as the attractor 

:= n n{M). (1) 

Note that (ZM = M \ Sing(A), and is thus not necessarily 
compact. 

Except in section [Vn| we shall generally omit reference to the 
singularity set Sing(A), but it is always understood that this 
set must be excluded when discussing the dynamics. 



A. Structural changes of strange attractors 

Figure [2] shows numerically-obtained attractors in phase 
spacefor pinbilliard dynamics in an equilateral triangle of side 
length 1, for several values of A. We see that trajectories 
from almost all initial conditions are observed to accumulate 
on fractal strange attractors. The main goal of this paper is 
the explanation of this phenomenology. The main structural 
changes observed are the following. 

For sufficiently small A, in fact A < ^, the attractor is a Can- 
tor set of lines, which "thickens" as A increases. For larger 
values of A, triangular- shaped gaps begin to appear in the at- 
tractor (b), which become larger as A increases. Afterwards, 
around A 0.68, the complete horizontal lines visible in the 
central bands of the attractor break up (c), which we refer to as 
band splitting. Soon thereafter, round A ^ 0.7, there is band 
merging to form just three bands (d). These bands shrink as 
A increases, concentrating around three angle values (e). Fi- 
nally, for A 0.98 and above, the attractor breaks into three 
distinct transitive regions (f). The enhanced online multime- 
dia version of this figure shows the evolution of the attractor 
for all values of A . These phenomena will be analyzed in turn 
in subsequent sections. 



B. Main theorem for A < ^ 

For A < ^, we are able to characterise the dynamics com- 
pletely via a symbolic model: 

Theorem 1 (Main Theorem). Let < A < ^. Then the at- 
tractor Ti of the pinball billiard map in the equilateral tri- 
angle is a transitive invariant set and is homeomorphic to 
[0,3] X C\Sing(A), where C is a Cantor set in R. More- 
over, there is an invariant ergodic measure on Tx, which is 
absolutely continuous with respect to the Lebesgue measure 
on horizontal lines. 
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FIG. 2. Phase space of pinbilliard dynamics in the equilateral triangle, depicting the main structural changes observed in the strange attractor 
as the contraction parameter A is varied. Each figure except (f) shows a single orbit of length 10^ starting from a random initial condition, after 
a transient of 10^ iterates were discarded to reach the asymptotic attractor. In (f), three such orbits are shown. (Enhanced online.) 



Section VII of this paper is devoted to a rigorous proof of the 
Main Theorem. 

Note that since all of these maps, restricted to the attractor, 
are conjugate to the same model, they are also conjugate to 
one another, i.e. the one-parameter family of pinbilliards is 
structurally stable for A < i . 



II. EQUILATERAL TRIANGLE 

For the pinbilliard map that we will study, we consider an 
equilateral triangle of side length 1 . The possible phase space 
coordinates are then M = (0, 3) x ( 



Uf^QM,-, where Mi := (/,/+!) x (- 



I , I) , which we split into 
I , I) is the part of phase 
space corresponding to the side /, with arc length in (/, /+ 1). 
(Note that the dynamics is not defined for phase space points 
whose arc-length parameter sits on a vertex, so these points 
are automatically excluded from the phase space.) When / 
denotes a side of the triangle, it must be read modulo 3, i.e. 
/ G Z3. In these coordinates, we will refer to vertical and hor- 
izontal intervals as being lines in phase space with constant 
first or second coordinate, respectively. 



A. Discontinuity curves 

Key objects in billiards with piecewise-smooth boundary are 
the discontinuity curves, which consist of those phase space 



points that are mapped to a vertex. For the equilateral trian- 
gle, the discontinuity curve consists of points {s^5{s)) G M 
with 



5{s) 



- arctan 



l+2(bJ 
v/3 



(2) 



where [^J denotes the integer part of ^ > 0. The function 5 is 
independent of A, and has image (— f , f )• As we shall see, 
this curve splits each Mi into two parts M^ and M^ , which are 
the domains of continuity of T^, see figure |3] 
Since the geometry of the equilateral triangle is sufficiently 
simple, we can compute explicitly the angular coordinate of 
7\ in terms of two affine maps: 



if 
if 



e>5{s), 
e<5{s), 



(3) 



where the functions (p±\ : R ^ R are defined by 

(p±i(e):=A(±f-e). (4) 
To see this, note that from a starting point (s^G) G M, on the 
side / := [^-J of the triangle, there are only two possible types 
of bounce, depending on 0: a right bounce, on side /+ 1, or 
a left bounce on side / — 1 ; these possibilities are separated 
by the point {s,5{s)) that hits the vertex. A straightforward 
computation verifies that the angular coordinate after the next 
collision is given by 0^(0), where the sign zG {+1,-1} de- 
notes a right or left bounce, according to the sign of the angle 
±1 between the sides before and after the bounce. 
The form of 7\ in ^ implies that locally, horizontal intervals 
in phase space are mapped onto other horizontal intervals with 
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FIG. 3. Schematic picture of the phase space of the equilateral tri- 
angle, showing the discontinuity lines (solid curved lines) and the 
resulting domains of continuity (shaded and empty). The Marko- 
vian region is shown, in addition to a rectangle (green and red open 
boxes) together with its images (green and red solid rectangles). The 
horizontal line shows the non-Markovian dynamics, with the open 
and closed circle indicating the orientation-reversing feature. The 
dashed vertical lines indicate the positions of the vertices. 



in general different angles, given by projections along a fixed 
direction from one side to another. 



B. Derivative 

To calculate expansion and contraction properties of the pin 
billiard map 7\ ' it is necessary to calculate its derivative at a 
point {so, Oq). The expression for this obtained for a general 
geometry^ reduces for polygonal billiard tables, which have 



only flat boundary components, to 



A B 
A 



with A -- 



cos 00 
COST]! ' 



B- 



COST]! 

(5) 

^0 denotes the distance between bounces in the billiard ta- 
ble and T]! the incoming angle at the next collision, as in 
secU 

The upper-triangular structure of the derivative of 7\ i^^" 
plies that horizontal tangent vectors of the form v = (a,0) 
are mapped back into horizontal tangent vectors with differ- 
ent base point. 



C. Periodic orbits 

Billiard maps do not have fixed points, and moreover in polyg- 
onal billiards neither are there period-2 orbits, unless two 
sides are parallel, in which case there is in fact a whole in- 
terval of such orbits. For pinbilliards with our rule, this still 
holds; thus, any periodic orbit in the equilateral triangle must 
have period at least 3. 



Lemma 2. Fix < A < 1. Then any periodic orbit ofT^ is 
hyperbolic of saddle type. 

Proof. Consider a periodic point (^o, ^o) with period n>?> and 
orbit {sj^Oj) := T^{so,6o) for j = l,...,f7. The periodicity 
implies that On = Oq. 

The derivative Z)(^q ^^^r^ is an upper- triangular matrix with 
eigenvector (1,0). The absolute value of the corresponding 
eigenvalue is 



cos Oq ^ cos 6i cos Oq cos a rji 

' cos T]„ COST],- COS rinfj[ COS T]/ 



(6) 



where r}i is the incoming angle at the ith collision. The left- 
hand side follows from the chain rule and eq. ^ after a rear- 
rangement; the equality follows from the pinbilliard collision 
rule Of = — At];, and using the fact that cos is an even function. 
Finally, we have cos0o/cosT]„ = cos(At]„)/cos(t]„), so that 
l^n = n?=i [cos(At]/)/cost];] > 1. Thus the line spanned by 
( 1 , 0) is the unstable eigenspace of the periodic orbit. 
Since the product of matrices of this form along the periodic 
orbit is also upper-triangular, its determinant is equal to the 
product of its diagonal entries and is also equal to the product 
of its eigenvalues. The other eigenvalue is thus A" < 1. Hence, 
there is one expanding and one contracting direction, and thus 
the periodic orbit is hyperbolic of saddle type. □ 
Note that this result holds for any periodic orbit of any polyg- 
onal pinball billiard, except the period-2 orbits mentioned 
above. 



D. Distinguished 3-periodic orbits 

For elastic billiard dynamics in an equilateral triangle, there 
are two distinguished period- 3 orbits, formed by the median 
triangle and an orientation. It turns out that we may continue 
these period-3 orbits for pinbilliard dynamics for all A < 1. 
They maintain their rotational symmetry and we have the fol- 
lowing lemma: 

Lemma 3. There are two hyperbolic period-3 orbits of saddle 
type for any A G (0, 1). Their local unstable manifolds are 
horizontal lines. 

Proof. The rotational symmetry of the distinguished period- 
3 orbit suggests that the outgoing angle must be a fixed 
point 6^ of the angular dynamics, so that (p^i{6^) = 6^. 
It is a straightforward computation that the unique solution 

is := e*(A) := f ^ = KA-A^ + A^ ). In fact, 

(p-i{—6^) = —6^is the counterpart with reversed orientation. 
A simple trigonometric calculation using the rotational sym- 
metry of the period-3 orbit shows that the point (^3 ( A ) , 0* ( A ) ) 
with 



*3(A):= 



\/3tan(0*)+3' 



(7) 



is periodic of period 3 for . 

The previous Lemma then shows that this period-3 orbit is 
always hyperbolic. 
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To see that the local unstable manifold of the period-3 point is 
a horizontal interval, let / C [0, 3] x {0h<} be a small horizontal 
interval that contains the period-3 point. Then (J) D /, so 
/ is a local unstable manifold of the period-3 point. □ 



III. MARKOVIAN REGION AND COMPLETE 
SYMBOLIC DYNAMICS 

We shall use the following feature of billiard dynamics in an 
equilateral triangle: Let Jq be a horizontal interval which cov- 
ers one side of the triangle: 

/e :=(/,/+!) x{e} CM,-, 

with / = 0, 1 or 2. There are two possibilities for the structure 
of the image T^{J), depending on the value of 6\ either it is 
a single horizontal interval or a union of two horizontal inter- 
vals, depending on whether the projection in the direction 6 
of the side / of the triangle does or does not hit a vertex; see 

fig-S 

When 1 1 > I , the image 7\ {J) is a horizontal interval con- 
tained in Mj, where j corresponds to the previous or next side 
in the table, depending on the sign of 6 ; and 7\ {J) is linearly 
contracted by a factor that depends only on 6 (and not on A), 
and which tends to 1 as 1 1 ^ | . 

On the other hand, when 1 1 < f the situation is very different. 
The interval / is now split by the discontinuity curve 5, and 
hence its image consists of two horizontal intervals, each of 
which completely covers its corresponding side j ± 1 . More- 
over, each is linearly expanded and maps across the whole 
horizontal extent of one of the other two sides. For this rea- 
son we call the subset [0,3] x (-f , f ) C M the "Markovian 
region". Note that these geometrical properties in the Marko- 
vian region in no way depend on A. 

The key observation is now that when A < ^ , the image (M) 
of the whole phase space is in fact contained in the Marko- 
vian region, and hence, roughly speaking, we can obtain a 
good partition in terms of future itineraries - i.e., whether the 
bounce is a left or right bounce at each iteration. This partition 



is described precisely in section VII A 



A. Angular dynamics 

The combinatorics of the horizontal intervals described in the 
introduction of this section depend on the dynamics of the 
angular coordinate, which we proceed to study in this section. 
The form of the maps 0± in ^ leads us to consider the fol- 
lowing geometric series with random signs. Let < A < 1 
and consider the following subset of real numbers: 

^(^) •= I f £ • ^ri e {1, -1}V^ > 1 I . (8) 



Suppose that we are given an angle 6 G C(A), that is 6 = 
f Lr=i^«^" for some choice of signs Zn G {+1,-1}. Then 
0±(e) gC(A). In fact. 



9±iW = 3 (±A-^^z.-iAM. (9) 



Thus the images <p±i(C(A)) C C(A) for both maps (p±i, and 
hence 

r;t([0,3]xC(A))c[0,3]xC(A), 

(except for points in Sing(A), where 7\ is not defined). Thus a 
good candidate for an invariant set is that it be contained in the 
set [0, 3] X C(A), which for A < ^ is a Cantor set of horizontal 
lines, and this corresponds with the pictures of the computed 
attractor shown in fig. [2ja). In section[VIl]we rigorously prove 
that the attractor is indeed precisely [0,3] x C(A) without the 
singularity set. 



B. Periodic orbits 

The geometric series that arise in the angular dynamics allow 
us to find periodic orbits with period m and particular itinerary. 
That is, suppose we are given a word w of length m > 3, given 
by w = {wj) with wj G {+1,-1} for j = 1, . . . ,m. Each wj 
specifies which kind of bounce, left or right, occurs at the jth 
iteration. If {s^^ 6^) is to be a periodic orbit of period m with 
this itinerary, then its angle must satisfy the equation 

= o • • • o (p^^ {0^) = - ^ (- 1)^- A^ + (-A)-e, 
^ i=l 

(10) 

where the second equality gives the explicit formula for the 
composition. Solving for 6^^ gives 



Ow — 



.7=1 



km 



k=0 



since the infinite series equals [1 — (—A)'"] ^ 

Setting bkm^j := (-1)^'"+^"" for any 7 = 1, . . . ,m and 

^ G N, we conclude that 

= ^ ^nA", 

and hence 6^ G C(A). 

When A < ^, all iterations of the interval [0,1] x {0^} are 
contained in the Markovian region. Hence we can partition 
this interval according to its future itineraries. Selecting the 
partition element corresponding to the itinerary w, we have 
T. 

itinerary (w, w, w, . . .) 
This is the basis of the argument for the rigorous proof of the 



^ (/vv) ^ Jw^ so that there is a periodic orbit with periodic 



We show in section VII that C(A) is a Cantor set for A < ^. 



Main Theorem presented in section VII Here we exhibit in 
fig.[5]two examples of periodic orbits found numerically using 
this method. 
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(a)Markovian (b)Non-Markovian 

FIG. 4. Images of horizontal intervals: (a) Markovian case, in which the other two complete sides are covered and there is linear expansion; 
(b) non-Markovian case, in which only one partial side is covered, and there is linear contraction. 



(a)Period 4 (b)Period 5 

FIG. 5. Periodic orbits (red) in the equilateral triangle for pinbilliard dynamics with A = 0.3: (a) period-4 orbit with itinerary 1,— 1) 

(b) period-5 orbit with itinerary The light grey lines show the complete regions with the same itinerary and angles. 



IV. NON-ACCESSIBLE REGIONS IN PHASE SPACE 

As discussed in the previous section, for A < ^ the attractor 
consists of complete horizontal lines, so that for a given angle 
6 G C(A), (almost) all positions in [0, 3] x {0} are attained in 
the attractor. For A > ^ , some of these horizontal lines are cut, 
giving rise to gaps, located close to the table vertices, that no 
longer form part of the attractor; see fig.|2jb). These regions 
are approximately triangular in shape, and change form as A 
varies. 

To explain this phenomenon, we say that a angle point (s^O) 
is non-accessible if {s^O) ^ T^{M)^ for some iterate n. For 

example, (0, 3) x |^ is such a non-accessible region for 

n = 1, along with its mirror image with negative angular coor- 
dinate. Obviously, such non-accessible regions cannot belong 
to the attractor. 

To avoid cumbersome notation, we will mainly work with in- 
tervals of angles and their dynamics under We denote 
by N := (— |, |) the whole range of available angles, and by 
K '= ( — f , f ) the angles corresponding to the Markovian re- 
gion. Furthermore, given an interval L = (^1,^2) and a scalar 
p, we denote pL := (p^i,p^2)- Recall that 0±i are affine 
maps of R that depend on A . 



Simple calculations show that the images 0+i {K) = (A | , A |) 
and 0-1 (^) = — 0+i(^) are disjoint intervals, both contained 
in ^N, and the gap between them is A^; that is, we have 
AA^ = 0+1 (K) U A/sT U 0-1 (K), together with two points at the 
boundaries of the intervals. 

If A < ^ , then we have that XN C K, and hence we can con- 
clude that the strip (0,3)xA^cMisa non-accessible region, 
since it is not contained in 7\ (M) . Iterating this procedure 
gives rise to the Cantor set C(A) discussed in the previous 
section. 

When A > ^, we now have the opposite case: AA^ D K, so 
there are points that may enter the middle gap. Since we have 
studied the images of points inside K under the action of 0+i 
and 0-1 , it remains to study the images of AA^ \ K. 
Let a := 0+i(+A|) |(2A - 3A^); by symmetry, we have 
0_i(-A|) = -a. For A > ^, we have a < A|, so that 
the image of AA^ \ K covers part of the previously-found gap. 
When a = 0, at A = I , the gap is completely covered by this 
image. 

The images of the interval \ K consist of two symmetric 
intervals, L := (a,A|) and — L. It is, however, necessary to 
study in more detail the position coordinate of the image under 
7\ of the strip (0,3) x L, since ?.N\K is no longer in the 
Markovian region, and so the image of horizontal intervals of 
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the form (/, / + 1) x do not completely cover a side. 
We thus only need to study the image of the rectangle R := 
[0, 1] X (f , A|) C MJ", where is in the domain of 
We have already shown that the image 7\ (^) ^ [1^2] x L. 
In order to obtain a better notion of the shape of the image, we 
have to be careful about the images of each of the four sides 
of the rectangle R. The map 7\ is, of course, not defined on 
vertical sides, but we can take its continuous extension. 
The bottom side of R, given by (0, 1) x {|}, is mapped onto 
the line (1,2) x {Af }, while the top side, (0,1) x is 
mapped onto a sub-interval (1,£) x {a}, reversing the orien- 
tation. The point £ = £{^) corresponds to the point where a 
ball hits the triangle's boundary if it starts on the vertex with 
angle A |. 

For the right side of R, it is necessary to consider vertical lines 
{1— e}x(|,A|)for small e > 0, since the image is not ac- 
tually defined for the right side of R. It is clear that the projec- 
tion on the angular coordinate of the image covers the whole 
interval (a, A |) , since the angles do not depend on the posi- 
tion s on the table. Furthermore, as e ^ 0, both ends of the 
curve converge to the vertex at ^ = 1 . So the continuous exten- 
sion of 7\ maps the right side of the rectangle onto the vertical 
interval {1} x (a,A|). 

Finally, the left side of R, given by {0} x (|, A|), is mapped 
to a curve joining the point (2,A|) and the point (^, a). 
In fact, this curve is given by Y{t) := (^(r),0+i(^)) for t G 
(f,Af), where 



1 + 73 tan(0 * 



(11) 



Note that only the endpoint ^ = ^(A|), and not the curve it- 
self, depends on A . 

To summarise, we have shown that 7\(^o^) is contained in 
Ml in the region bounded below by the curves max{a, Y{t)}, 
and above by the line A | . 

The other way that trajectories can be injected into Mi is if 
they come from M2 . By symmetry, the relevant image of M2 
is bounded by the curves — A| and min{ — a, 7(^)}. This is 
shown in figure [6j and explains the presence of the triangular- 
shaped gaps in the attractor. 

As we have seen, when A > | , the angular parts of the images 
overlap. In fig. [6|c), we see that the complete images in fact 
touch. Although this makes possible that the central bands 
of the attractor join as soon as A > |, in fact we observe nu- 
merically that this occurs later, around A :^ 0.7. This can be 
understood by tracing the evolution of the non-accessible re- 
gions in time and finding the limits of the region around 0=0 
which is correspondingly excluded from the attractor. 



V. ISOLATION OF PERIOD-3 ORBIT FROM THE 
ATTRACTOR 



Figure [7] shows the period-3 orbit together with the attractor. 
From the figure, we see that the period-3 orbit is embedded in 
the attractor before the bands split, and in particular that its lo- 
cal unstable manifold is a complete horizontal line contained 
in the center of the band. After the band splits, however, the 
period-3 orbit becomes isolated from the strange attractor, and 
there is no longer such a central line. 

This phenomenon is due to a homoclinic bifurcation: for 
A = 0.55, the stable and unstable manifolds of the period- 
3 meet, forming a homoclinic point, as shown in fig. [Sja). 
Thus any neighbourhood of the period-3 orbit contains recur- 
rent points, and hence the period-3 point is not isolated from 
the dynamics. 

By A 0.75 the stable and unstable manifolds no longer meet, 
and this recurrence is destroyed, allowing the period-3 orbit to 
become isolated. 

To obtain these figures, the unstable manifold is obtained by 
iterating a local horizontal interval around the period-3 point. 
The stable manifold, on the other hand, is approximated us- 
ing the following escape-time algorithm. For a given point 
in phase space to belong to the local stable manifold of the 
period-3 orbit, its future iterates must follow the same dynam- 
ics as the period-3 orbit, e.g. only have right bounces. We 
thus iterate each point in a grid on the domain a certain num- 
ber n of times; if the point has only right bounces up to this 
time, then the initial condition is marked as forming part of 
the approximate local stable manifold of order n. Of course, 
this does not give the precise shape of the stable manifold, but 
rather a region in which it must be contained. In particular, 
this guarantees that there is no homoclinic point if the unsta- 
ble manifold does not cross this region. 
This phenomenon also provides evidence of the relatioJ^ 
(Section III. 3) between expansion in the unstable direction 
and fractioning due to the set Sing(A). The expansion in the 
unstable direction of the period-3 orbit is given by [^^^y]^, 
which decreases as A increases. 



VI. NON-TRANSITIVE BEHAVIOUR CLOSE TO A = 1 

For A 0.98, we observe that there are three distinct transitive 
parts of the attractor, i.e. such that it is impossible to jump 
from one to the other; see fig. [9] 

The reason for this phenomenon is the following. After 
enough iterations, the system arrives close to a vertex, with 
an angle close to | , so that it collides with the opposite side at 
an almost perpendicular angle, and is reinjected back into the 
same corner. This acts as a trapping region which prevents the 
system from escaping to one of the other parts of the attractor. 



In 10 the basins of attraction of each of the three symmetri- 
cal copies of this region are shown, indicating a fractal basin 
boundary structure. 



Recall from figure |2jb) that a splitting of the central bands of 
the attractor is observed around A 0.68. This is related to 
the local stable and unstable manifolds of the distinguished 



period-3 orbit discussed in section IIP 



VII. PROOF OF MAIN THEOREM 

In this section we give a rigorous proof of the Main Theo- 
rem. This will be done in two parts, treating separately the 
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(a)A = 0.3 (b)A = 0.5 (c)A = 0.75 



FIG. 6. Action of on AMq for different values of A. The Markovian region is shaded in green and the non-Markovian region in blue. Darker 
shading is used for the images, shown for 5 G (1,2); lighter shading is used for the domain, shown for 5 G (0, 1). The attractor is shown in red 
dots. 
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(a)A = 0.5 (b)A = 0.68 (c)A = 0.8 

FIG. 7. Phase space of pinbilliard dynamics in the equilateral triangle, with the period-3 orbits shown as crosses. 



topological and the ergodic properties, although both are con- 
sequences of the existence of a topological conjugacy of the 
pinball billiard map in the attractor with a symbolic dynamics 
model. 

Let us recall that the attractor of 7\ (or maximal forward 
invariant set in M) is defined by •= ^n>oTx{M) C M. 
Roughly speaking, the symbolic model will be a shift- like 
map in the space of bi-infinite sequences of ±1, but inter- 
changing signs in a certain part of the sequence at each it- 



eration; this map is defined in eq. p4| ). This map possesses a 
natural ergodic measure which will be pulled down to an in- 
variant, ergodic measure for the pinball billiard map using the 
above-mentioned conjugation. 

The Markovian property that guarantees uniform expansion 
along certain horizontal intervals is the cornerstone in the con- 



struction of the conjugation, discussed in sec. Ill and con- 
structed in section [VII A| below. In some cases, for instance 
for < A < ^ , the attractor is contained in the region where 
this property holds, allowing us to prove the existence of a 
dense orbit. 

It is interesting that the set of geometric series with random 
signs is intimately related with the action of the collision rule 
of the pinball billiard map on the angular coordinate. 
Let < A < 1 and recall the set C(A) defined in ([s]): 

C{X) := 1 1 £ z„A" : z« e {1,-1}V« > 1 1 . 



This is a subset of the interval | [— , jzj]- It is a Cantor set 
if A < ^, and all numbers in C(A) then have a unique expres- 
sion in terms of a signed geometric series, whereas for A > ^ 
it is a full interval, which grows to fill up R as A ^ 1. 
To prove the Cantor structure of this set, note that if A < ^ 
then for each integer m > 1 we have 



j=m+l 



71 1- 

3 T 



-2A 



A 



A^>0. 



Hence, there are no points of C(A) in any open interval of the 
form {s — 5m,s -\- 5fn), for any partial sum s = j Y!J=i ±>^^- 
In the same fashion, the length of each bounded connected 
component of the complement of the union of these intervals 
is of order A^+^. 

On the other hand, if A < | then C(A) C [-f,f]. Re- 
call that the set C(A) is invariant under both maps (p±i, i.e. 
(p±i(C(A))cC(A). 

To observe the relation between the geometric series and the 
pinbilliard, consider the topological space 

I.t:={{zn):zne{h-l}\/n>l}, 

with the standard product topology, and the modified shift 
transformation — cr : ^ I^J given by —a{{zn))'= (— z^+i)- 
In these terms, if < A < ^, then both maps (p±i are branches 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



(a)A = 0.55 (b)A = 0.75 

FIG. 8. The region Mq with non-negative angles, showing the strange attractors (blue points), one point in the period-3 orbit (crossed circle), 
several iterations of the unstable manifold of the periodic point (thick black line), and a region in which its local stable manifold must be 
contained (red area) for (a) A = 0.55 and (b) A = 0.75. In (a), there is a homoclinic point (intersection of the stable and unstable manifolds), 
marked with a small red circle. A number of iterations n = 30 was used. 




FIG. 9. Three non-communicating transitive parts of the attractor, shown in different colors for A = 0.98, in (a) configuration space and (b) 
phase space. 



of the inverse (— cr(z)), where the map E;^ : 
is defined by 



■C{X) 



(12) 



n=\ 



The condition < A < ^ guarantees that E^^ : C(A) is 
well-defined. 

On the other hand, to any point in phase space having an entire 
well-defined future orbit there corresponds a sequence (w„) G 
, such that w„ = +1 if the ni\\ iterate is a right bounce, or 
w„ = — 1 if it is a left bounce. So, if {sq^ Qq) G M\ Sing(A), 
then for any ^ > 1 we have T^{so^ Gq) = {sn^On), where 



that the tail of the itinerary is either Wk^j = (— 1)^ or (— 1)^+^ 
for J > 1. 

Recall that for any < A < 1 we have Tx{M) c (0,3) x 
(-Af ,Af), and note that if < A < ^ then C(A) C 
(— A|, A|). For such A, the attractor inherits the geome- 
try of C(A), as shown in the following lemma. 

Lemma 4. 7/0 < A < ^ then c (0,3) x C(A). 

Proof. Consider any point (q^O) G r;^- Then there exists a 
point (pojOo) ^ M\Sing(A), and an increasing sequence of 
integers rik such that (;?„^, := T^^{po, Oq) ^ {q,0) 
as ^ ^ oo. Using the expression in ([13]), we have, for any 
k> 1, that 



e. = 9w.o---o(p^^(eo) = -X(-i)^-V,_,-+iA^+(-A)"eo. 

(13) 

In this way, we can define the function : ^ 
J^{P) := {wn) e Zj for P er^. if a point hits a vertex of 
the table at some iterate ^ > 1 then one cannot compute its 
itinerary after the ^th step. In this case, we take the convention 



n -- 



£(-l)^-V„,_,+iA^- 




for some positive constant C and any selection of ±1 in the 
series on the right-hand side. Thus Q G C(A). □ 
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FIG. 10. Basins of attraction of the three transitive parts of the attractor, shown in three different colors, exhibiting a fractal basin boundary, 
(a) Region around one part of one of the transitive parts of the attractor, centered at (0.5, 0). (b) Detail of the basin boundary. 



A. Region with a Markovian property 

The discontinuity curve 5, defined in ([2]), does not depend 
on the parameter A, and its image is contained in the rect- 
angle (0,3) X (— f,f). Furthermore, phase space is natu- 
rally partitioned into horizontal intervals of the form /(/, 0) := 
(/,/+!) X {0} for / e Z3 and e (-f , f ). We denote this 
partition by 

We define ^0 C ^ as ^0 := {J{i.O) e / : \Q\ < f }. 
The region of phase space covered by i-e. the union of 
all /(/, 0) G ^ such that \Q\ < |, has a non-uniform ex- 
panding behavior along horizontal lines. To see this, note 
that any / G ^0 is partitioned into two non-trivial subinter- 
vals / = /(/,e) = Vq[JR'q, where Uq := {i,de] x {0} and 
R^Q := 1) X {0}, and is the first coordinate of the 

intersection point of the image of the curve 5 with the interval 
J{i,e). So, the ^Qiy = {Uq,R'q :/GZ3,|e| < f} is a refine- 
ment of such that any atom of y is almost contained in 
some continuity domain of 7\ • Precisely, given any S ^ ^ , if 
we denote by int(5') the interval S without its extreme points, 
then we have that int(5') is contained in some continuity region 

The Markovian property of the partition is the following: 
for any S G the restriction of 7\ |int(^) extends to a contin- 
uous and surjective map fs : S ^ J, for some J e In fact, 
each of these maps is a linear projection between sides of the 
triangle, along some fixed direction, and hence affine. Fur- 
thermore, for any 5" G ^ the derivative is constant and can be 
computed as {fsY = l^"!"^ > 1, where l^"! denotes the length 
of the interval S. Note that {fj{i^Q))' ^ 1 as ^ |. These 
properties do not depend on A . 

If we fix some < A < ^ , then the set is properly contained 
in the Markovian region. That is, there is e > such that 

r;iC(0,3)x(-f + e,f-e)cUye^o-^' 
where equality is excluded. 

So for any 5" G ^ we have that fs{S) G ^q. Moreover, there 
is a uniform expansion along the horizontal direction: in fact, 
there is p = p(A) > 1 such that {fs)' > p, for any S C /(/, 6) 
such that |e| < Af . 

These properties allow us to describe more precisely the dy- 
namics of the pinbilliard map in terms of symbolic dynamics. 



Let us first state the following lemma: 

Lemma 5. Fix some < A < ^. For each (/, 0) G Z3 x C(A) 
there is a continuous map P{i^ 0) : ^ /(/, 0) G ^0 such 
that y oP = ld. 

Proof. Fix some < A < ^, and consider any(/, 0) G Z3 x 
C(A). Let Jo = J{i,0) G ^q. In order to define the map 
P(/, 6) : ^ Jq, consider a point w = (w„) G . We take 
the first k letters of w and will find a closed interval /4 = 
//(wi , . . . , Wy^) cJq such that, for any Q e H^, the first n coor- 
dinates of y{Q) G Zj are precisely wi , . . . , w^. If HkH^ C Jq 
is a unique point, then the map will then be well-defined. 
To define these intervals we shall proceed inductively on k: 
Recall that Jq G and so it splits into two subintervals, 
Jo=LUR, for some L,Rey\ then let Hi=Rifwi = 1, or 
H\ = L if wi = — L It is easy to see that the first coordinate of 
J^(2) is wi , for any 2 G if 1 . Notice that fj^ (if 1 ) =: /i G ^o- 
Now let us assume that H^-i C /q is defined and satisfies 

4:=/4_iO---o//iO//o(^^-i) ^ /o. 

Then //. = L^URk for some L/., R^ G and so, in the same 
way, we can define = {fj^_i o • • • o//^)"^ (A), where either 
A = Ly^, or A = i?yt, depending on w^. 

We have that |i^^| ^ as ^ ^ 00. In fact, all Jk=J{-, %) G 
^0 involved in the construction are such that | % | < A | < | , 
and hence {fjj > p > L So \Hk\ < p~^, for any ^ > 1, and 
therefore, there is a unique P = P{i^ 0){w) G n^>i and we 
have y{P) = w. 

This map is continuous by construction and it is injective at 
all points, except those with a periodic tail (see the definition 
of the map J^), which correspond to points that eventually hit 
some vertex of the table, which are mapped to the same points 

in/o- n 
As a consequence of the construction in the proof of the pre- 
vious lemma one obtains the following: on any /:=/(/, 0) G 
^0 such that JPiTx 7^ 0, that is G C(A), any itinerary 
in Zj is followed by a unique point in /, with exception of 
those itineraries with periodic tail, which correspond to the 
same orbit that eventually hits a vertex. In fact, such excep- 
tions correspond to points in the discontinuity set Sing(A). 
Hence, / C ' and this proves that is homeomorphic to 
(0,3)xC(A)\Sing(A). 
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Thus, we can gather together the maps obtained in Lemma |5] 
to create the following transformation that depends only on A : 

/\ :Z+xz+xZ3^r;t; 

/\(w,z,/):=P(/,^aW)W. 

Observe that this map is continuous, and if (p^O) G (0,3) x 
C(A) then the inverse of is defined by P^^{p^6) = 
{y{p^6)^E^^{6)^[p\), where [p\ denotes the integer part of 
p. Hence, is a homeomorphism in P^^ {T^ \ Sing(A)) and 
conjugates the pinball billiard map to the following symbolic 
model: 



and iterate k-\-2 times to obtain the claim. In fact. 



[(w,z),/] ^ [((t(w), (wi, -Zi, -Z2, . . .)),i^wi] 



(14) 



where O" : ^ ^2 standard shift, and the addition op- 

eration in the last coordinate is in Z3 . 

This complete description of the dynamics of the pinball map 
allows us to prove transitivity in ' we shall see in the next 
Lemma. 

Lemma 6. The pinbilliard map Ti is transitive. 

Proof. Take any two non-empty open subsets U and V of • 
We shall prove that there is some > such that {U){~^V ^ 
0. First note that there is some interval J\ := J{i\^B\) G ^0 
such that /i n is a non-trivial interval, since U is open. Now, 
the construction in the proof of Lemma [5] allows us to say 
that there is some interval H\ C JiDU and Ni >0 such that 
(Hi) =J2e J^. Therefore, there are (/2, ^2) ^ Z3 x C(A) 
such that Ji= J(ii^Qi). 

On the other hand, there is some interval J{i^jO) G such 
that /(4, 0) ny 7^ for some (/*,§) G Z3 x C(A), since V is 
open. 

Now we claim that for any £ > there is some N2 > and 
an interval H2 C J2 such that T^^ {H2) = ^3 G and hence 
^3 = J{i3^^3) for some {is ^63) G Z3 x C(A), and such that 
/3 = and 1 03 — 1 < e. If this claim holds and we take £>0 
small enough, it implies that T^^^^^{Hi) flV 7^ 0, proving 
transitivity. 

To get the claim denote (ai , <22, . . .) = ^^H^) ^ ^2 • some 
^ > 1 such that < £. We will select a point w G 
specifying just the first k coordinates in the following way: let 



w:={{-lf-'at,{-ir-'at-i 



-a2,ai,wyt+i,...) G Zj. 



If wQlQtz = £^^{62), then by the formula in dl4l one gets 



a+2 



(w,Z,/2) = (^(T^^^(w), (ai,...,a;^,Z?i,Z?2,Zl,...) 



h^Y.^{-iy-'as-b2^hj. 

If these numbers differ only by 1, then the same argument 
applies by adjoining only one /?i G {±1}. □ 
This completes the proof of the topological part of Theorem 

m 



B. Invariant measure on 

The space x x Z3 has a natural probability measure 
/i, defined on Borel sets, which is invariant under the map 
in ( p4| ). This measure is in fact a Bernoulli measure, and is a 
product /i = V X V, where v is the invariant Bernoulli measure 
for the standard shift a : Zj ^ Zj , weighted in order that 
;U(Z+ X Z+ X {/}) = ^, for / G Z3. 

On the other hand, one can define a measure fl^ on Borel sets 
of R using the maps E^^for < ^ < 1, in the following way: 
given any measurable set A c R, let 

Ma(A):=v(£^-1(A)). 

Properties of the measure fl^, in particular for which val- 
ues it is absolutely continuous with respect to Lebesgue 
measure, have been studied since 1920 by Khintchine and 
Kolmogorofl^, and several interesting results have been de- 
veloped by Erdos, Garsia, Pollicott and Solomyak during the 
20th century; se d^^f^^i In particular, Jensen and Wintner irt^ 
proved that the measure jl;^ is either purely singular or ab- 
solutely continuous with respect to Lebesgue measure in the 



interval | 



1-A ' 1-A 



IfA< 



1 



2 then the measure jlp^ is purely 

singular, since it is supported on the Cantor set C(A). 
So, if < A < ^, the measure fl gives us an invariant measure 
for 7\ on • Let m be a measure on Borel sets of A C M 
given by 

m{A) :=Ai(/\-'(Anr;i)). 

Theorem 7. If{)<X < ^, then the measure m is an invariant 
ergodic measure for the pinball billiard map Ti. Moreover, 
the conditional measures m\j, for J G ^0 such that / fl 7^ ®' 
are absolutely continuous with respect to the Lebesgue mea- 
sure on the interval. 



T^(w,z,i2) - (<7^(>^), {ai,...,ak,z\,...), ^2 + ^( 1)^ ^cis). Proof By construction, the measure m is 7^ -invariant, since it 



s=\ 

Observe that \Q - Ex{ai,. . . ,ak,z\,. . .)\ < < £. So, 
in order to obtain the claim it is left to prove that i2 + 
Y!l=\{—^y~^ = /*niod3. However, this is not always true. 
For instance, let us assume that they differ by 2. Then, to get 
the claim we may need to change the itinerary w, by adjoining 
two signs b\ and /72 ^ {±1} at the beginning: 



w - 



{b2,bi,{-lf ^ak,{-lf ^ak-u...,-a2,ai,...) 



\k-2 



is topologically conjugate to the map in ([14]), and hence it in- 
herits all its ergodic properties. On the other hand, conditional 
measure on horizontal lines are absolutely continuous with re- 
spect to the Lebesgue measure, since the map 7\ restricted to 
/ G ^0 that Jnr^ / is uniformly expanding with a factor 
of p > 1. Standard methods (Theorem 2.1) complete the 
proof. □ 
Theorem [7] proves the ergodic part of Theorem [T] and we are 
done. □ 
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VIII. CONCLUSIONS 

We have built up a detailed picture of the dynamics of non- 
elastic pinbilliards in an equilateral triangle as the contraction 
parameter A varies between and 1. Away from these lim- 
its, the pinbilliard has dominated splitting, and we have es- 
tablished, rigorously for A < ^ , an intricate limiting dynam- 
ics consisting of chaotic strange attractors which evolve as A 
varies, and which undergo qualitative changes at certain well- 
defined parameter values. 

Preliminary results indicate that non-equilateral triangle ta- 
bles exhibit similar phenomena, except that the limiting elas- 
tic dynamics is then expected to be generically mixing, which 
seems to be reflected in the attractors found when A close to 
1, which fill out more and more of phase space as A ^ 1, in- 
stead of the shrinking and non-transitive behaviour we have 
discussed here. Other polygonal billiard tables also seem to 
exhibit similar behaviour, with additionally the possibility of 
families of attracting period-2 orbits if the table has two par- 



allel sides. 
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